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Abstract 

In this presentation we overview some recent results on biased tracer diffusion in lattice 
gases. We consider both models in which the gas particles density is explicitly conserved 
and situations in which the lattice gas particles undergo continuous exchanges with a 
reservoir, which case is appropriate, e.g., to adsorbed monolayers in contact with the 
vapor phase. For all these models we determine, in some cases exactly and in other ones 
- using a certain decoupling approximation, the mean displacement of a tracer particle 
(TP) driven by a constant external force in a dynamical background formed by the lattice 
gas particles whose transition rates are symmetric. Evaluating the TP mean displacement 
explicitly we are able to define the TP mobility, which allows us to demonstrate that the 
Einstein relation between the TP mobility and the diffusivity generally holds, despite the 
fact that in some cases diffusion is anomalous. For models treated within the framework 
of the decoupling approximation, our analytical results are confirmed by Monte Carlo 
simulations. Perturbance of the lattice gas particles distribution due to the presence of a 
biased TP and the form of the particle density profiles are also discussed. 



1 Introduction. 



Transport in media with quenched or time-dependent disorder has received much interest due 
to the considerable fundamental and technological importance of the problem [1-8]. Transport 
processes underlay a variety of chromatographic separation techniques [4-6], extraction of 
oil from porous rocks [4-6], excitations quenching in amorphous solids and doped crystals 
[6, 7, 9-12], conductivity of complex media and permeability of disordered membranes [8, 13]. 
Molecular diffusion measurements in disordered media often serve as a tool in characterizing 
the internal geometry over a range of molecular and macroscopic length scales [4-7]. In 
addition, transport processes control the kinetics of chemical reactions in disordered media 



Following motion of individual particles in disordered media one most often encounters 
behavior characterized by a mean-square displacement < R 2 (t) > growing as 



i.e. linearly in time with D being the "diffusion" coefficient, which depends on the peculiar 
features of the disordered medium. Such type of dynamics is ubiquitous due to its inherent 
relationship to the central-limit theorem, which makes it a paradigm of stochastic processes. 

One the other hand, there was a growing evidence for other laws of random transport, 
both slower and enhanced relative to the time dependence in Eq.(l). Generally, one finds 



Case with a < 1, known as the dispersive transport or subdiffusive regime, has been in the 
focus of many theoretical and experimental studies (see, e.g., Refs.[l-4, 12]) and has been 
attributed, for instance, to random walks on fractal structures [1-4, 12] or to diffusion in 
presence of temporal traps [1,12,14]. A few stray examples of such a transport are time 
of flight measurements in amorphous semiconductors [14,15], defect diffusion in relaxation 
studies [14], excitation migration on the disordered array of donor centers (see, e.g., Ref.[16]), 
tracer diffusion in a linear array of convective rolls [17] or in low-dimensional lattice gases 
[18-23], or the so-called Sinai diffusion with extremely confined particle trajectories such that 
< R 2 (t) >~ ln 4 (t) [24-26]. Enhanced diffusion regime, a > 1, is observed in turbulent 
diffusion (a = 3), as well as in some other instances (see, e.g., Ref.[27]) for some review). 

Suppose now that in the disordered medium in question, in which, in absence of any 
fields acting on the species whose motion we follow (say, tracer particles), we find that their 



[9-12]. 



< K 2 {t) >~ Dt, 



(1) 



< R 2 (t) >~ t a , with a / 1. 



(2) 
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mean displacement is zero, while the mean-square displacement < R 2 (i) > obeys Eq.(l) or 
the anomalous diffusion law in Eq.(2). Suppose next that we manage to bias somehow the 
particle motion, say, by charging the tracer particles only and putting the whole system in 
a constant electric field E. This situation is, of course, quite interesting in its own right 
and is appropriate to charge transfer or "dynamic directed percolation" in time-dependent 
inhomogeneous medium [28-31], such as, e.g., certain biomembranes [32], solid protonic con- 
ductors [33], oil-continuous microemulsions [34-37] or polymer electrolytes [38-40]. In such 
a situation, the tracer particles will move preferentially in the direction of the field and will 
attain a non-zero mean displacement and a non-zero mean velocity V(t). One may define 
then the time-dependent tracer particle mobility as 

H(t) = lim ^1, (3) 

while the " diffusivity" D(t) can be extracted from the behavior in absence of the field (Eqs.(l) 
or (2)) via: 

where d is the spatial dimension of the system. 

Now, the question is whether the mobility fj,(t), calculated from the tracer particle mean 
displacement in the presence of an external electric field, and the diffusivity D (t), Eq.(4), 
deduced from the tracer particle mean-square displacement in the absence of the field, obey 
the generalized Einstein relation of the form 

u(t) = (3D(t), (5) 

where (3 denotes the reciprocal temperature? 

The answer is trivially positive, of course, for isolated non-interacting particles performing 
conventional Brownian motion. For particle diffusion in disordered medium the validity of the 
relation in Eq.(5) is much less evident. On the one hand, it has been found that Eq.(5) holds 
for the tracer particle motion in a one-dimensional (ID) gas of hard spheres, whose whole 
dynamics consists of pairs of neighboring particles interchanging velocities at each collision 
and the initial velocity distribution being an equilibrium one [41]. It holds also for the tracer 
particle diffusion in a ID hard-core gas with diffusional dynamics when the ID lattice is a ring 
of a finite length [42], as well as for infinite ID lattices with non-conserved [43] and conserved 
particles number [44-46]. Remarkably, in the latter case Eq.(5) holds for t sufficiently large 
despite the fact that both u{t) and D (t) do not attain constant values as t — > oo but rather 
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are decreasing functions of time [44-46]. In addition, the validity of the Einstein relation has 
been corroborated for the charge carriers in semiconductors [47] and for polymeric systems 
in the subdiffusive regime [48, 49]. 

On the other hand, it is well known that the Einstein relation is violated in some physical 
situations; for instance, it is not fulfilled for Sinai diffusion [24-26] describing Brownian motion 
of a test particle in presence of a quenched random force, and for the Scher-Lax-Montroll 
model of anomalous random walk [50] (see also Refs.[51] and [52] for some other examples). 
Apparently, the Einstein relation is also violated for random walk in disordered lattices at their 
percolation threshold - the celebrated model of the "ant in the labyrinth" of de Gennes [53]. 
Here, in absence of external force acting on the random walker, its mean-square displacement 
< X 2 (t) > obeys Eq.(2) with a < 1 [53, 54]. In presence of external field, which produces a 
bias in the random walk making the walker more likely to step along the field than against it, 
the behavior appears to be rather complex. As a matter of fact, the bias has a dual effect - it 
induces drift in the direction of the field, but also creates temporal traps, such as in dead-end 
branches, to escape from which the walker has to move against the field. One then finds that 
the drift velocity is not a monotonic function of the bias [55, 56] and moreover, it has been 
argued that the drift velocity vanishes once the bias exceeds a certain threshold value [56] 
(see also [57]). These theoretical arguments have been supported by exact calculation of the 
mean drift velocity on a random comb and randomly diluted Bethe lattice [58] (see, however, 
Ref.[59] for some objections). Therefore, a small value of external bias gives rise to a mean 
displacement in the direction of the field that increases linearly with time, and the mean 
velocity tends to a constant. This asymptotic value of the mean velocity is proportional to 
the field for small fields [57], which implies that the mobility tends to a constant as t — > oo. 
Contrary to that, since in the absence of the field the walker motion is subdiffusive [53, 54], 
one has that here the diffusivity D(t) should vanish as t — > oo. Therefore, if the results of 
Refs. [55-58] are correct, they would imply that the Einstein relation is not valid for diffusion 
in random lattices at percolation threshold. 

Hence, in principle, it is not a priori clear at all whether Eq.(5) should be valid in any 
case. 

In this presentation we overview some recent results on the generalized Einstein relation 
in Eq.(5) for tracer particle (TP) diffusion in dynamical disordered environments, as exem- 
plified here by non-interacting lattice gases of hard-core particles. The paper is structured as 
follows: In Section 2 we introduce the model, basic notations, write down dynamic equations 
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describing the time evolution of the system and discuss the way of their solution. In Section 3 
we first recall the properties of an unbiased TP diffusion in a ID hard-core lattice gas. Next, 
we turn to the case of a biased TP in an inert lattice gas and present explicit results describing 
the TP mean displacement and the mobility. We demonstrate then that the Einstein relation 
in Eq.(5) holds exactly despite the fact that in this case diffusion is anomalous. Further on, 
in Section 4 we consider the case when the initial particle distribution on a ID lattice is inho- 
mogeneous and characterized by an S*-shape, "shock" -like density profile. We find that when 
the TP is attracted towards the particle phase with higher density, at a certain critical value 
of the constant attraction force the TP will not move, on average. Assuming the validity 
of the Einstein relation, we find the second moment of the TP displacement. This heuristic 
analytical result is confirmed numerically which signifies that the Einstein relation holds ap- 
parently for such an inhomogeneous situation. Next, in Sections 5 and 6 we analyse the forms 
of the biased TP terminal velocity in one-, two- and three-dimensional hard-core lattice gases 
undergoing continuous particles exchanges with a reservoir. We determine the TP mobility 
and assuming the validity of the Einstein relation, obtain the TP diffusion coefficient D in the 
unbiased case. Results for ID systems are confirmed by numerical Monte Carlo simulations. 
We also show that in two- and three-dimensional systems in the conserved particle density 
limit (when particle exchanges with the reservoir are forbidden) our predictions for D coin- 
cide with classical results of Nakazato and Kitahara [60], which are known to be in a very 
good agreement with the numerical data [61]. In Section 7 we present an exact solution for 
the biased TP mean displacement in a two-dimensional (2D) lattice gas containing a single 
unoccupied site. On comparing our results with earlier results by Brummelhuis and Hilhorst 
[62] obtained for the unbiased TP dynamics, we infer that also in this case, in which the 
TP diffusion is anomalously, logarithmically confined, the Einstein relation does hold exactly. 
Finally, we conclude in Section 8 with a brief summary of results and discussion. 



2 The model and basic equations. 

Consider a d-dimensional regular hypercubic lattice of spacing a each site of which is brought 
in contact with a reservoir containing identic, electrically neutral particles - a vapor phase 
(Fig.l), maintained at a constant pressure. We suppose next that the reservoir particles 
may be created (adsorbed) on any vacant site at a fixed rate f/r*, which rate depends 
on the reservoir pressure and the energy gain due to the "creation" event. Further on, 
the particles may move randomly along the lattice by hopping at a rate l/2dr* to any of 



4 



2d neighboring lattice sites, which process is constrained by hard-core exclusion preventing 
multiple occupancy of any of the sites. Lastly, the particles may spontaneously disappear 
(desorb) from the lattice at rate g/r*, which is dependent on the barrier against desorption. 
Both / and g are site and environment independent. 




Figure 1: Two-dimensional realization of the model: Adsorbed monolayer in contact with a vapor 
phase. Gray spheres denote the monolayer (vapor) particles; the smaller black sphere stands for the 
driven tracer particle. 

To describe the occupancy of lattice sites, we introduce a time-dependent variable ry(R), 
which may assume two values: 



r?(R) 



1, if the site R is occupied, 
0, if the site R is empty. 

Note that the local variable r/(R) can change its value due to creation of particles (adsorbtion), 
desorption and random hopping events. Note also that random hopping events do conserve 
the total number of adsorbed particles and hence, the average particle density p s (t). On 
the other hands, creation/desorption processes do change ??(R) locally such that the total 
number of particles on the lattice is not explicitly conserved when such two processes are 
allowed. However, the mean particle density, p s (t) =< rj(R) >, approaches as t — > oo a 
constant value 

This relation is well known and represents the customary Langmuir adsorption isotherm [63]. 
We finally remark that in the analysis of the stationary-state behavior, we can always turn 
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to the conserved particles number limit by setting / and g equal to zero and keeping their 
ratio fixed, i.e. supposing that f/g = p s /(l — p s )- This limit will correspond to the model of 
biased tracer diffusion in a hard-core lattice gas with fixed particles density p s , and will allow 
us to check our analytical predictions against some already known results [60-62, 64, 65]. 

Now, we introduce a bias into the system. Contrary to the well-studied situation in which 
all particles are subject to a constant force (see, e.g., Rcfs. [66-69] and references therein), we 
will focus here on essentially less studied case when only one of particles is subject to external 
bias while the rest have symmetric hopping probabilities. This particle may be thought of as 
a probe designated to measure the frictional properties of symmetric non-inteacting lattice 
gas. To this end, at t = we introduce at the lattice origin an extra hard-core particle, 
whose motion we would like to follow and which will be called in what follows as the tracer 
particle - the TP; position of this particle at time t is denoted as Rj r . We stipulate that the 
TP is different from other particles in two aspects: a) it can not disappear from the lattice 
and b) it is subject to some external driving force, which favors its jumps into a preferential 
direction. Physically, (a) be realized, for instance, if this only particle is charged and the 
system is subject to a uniform electric field E. We suppose here, for simplicity of exposition, 
that the external force E is oriented according to the unit vector ei. 

The dynamics of the biased TP is defined in the usual fashion [23]: We suppose that the 
TP, which occupies the site Rt r at time t, waits an exponentially distributed time with mean 
r, and then attempts to hop onto one of 2d neighboring sites, R tr + e^, where e u are 2d unit 
vectors of the hypercubic lattice. In what follows we adopt the notation v = {±1, ±2 . . . , ±d}, 
where ei will denote the direction of the external force E. Next, the jump direction is chosen 
according to the probablity p„, which obeys: 



exp 


f(E-e^ 




exp 





(7) 



where j3 is the reciprocal temperature, (E • e) stands for the scalar product, the charge of 
the TP is set equal to unity and the sum with the subscript \i denotes summation over all 
possible orientations of the vector e^; that is, fi = {±1, ±2 . . . , ±d}. 

After the jump direction is chosen, the TP attempts to hop onto the target site. The hop 
is instantaneously fulfilled if the target site is vacant at this moment of time; otherwise, i.e., if 
the target site is occupied by any adsorbed particle, the jump is rejected and the TP remains 
at its position. 

To close this subsection, several comments on the relevance of the model to some physical 
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situations are in order. As we have already remarked, in the conserved particles number limit 
the model under study reduces to the model of biased tracer diffusion in hard-core lattice gases, 
which is related to a wide variety of physical situations, including, for example, diffusion of 
interstitials in crystals. A detailed review can be found in Ref.[61]. In two-dimensions, in 
the non-conserved particles number case it describes tracer diffusion in adsorbed monolayers 
in contact with the vapor phase [70] and allows to determine, for instance, the analog of the 
Stokes formula for such monolayers and consequently, to elucidate their intrinsic frictional 
properties. On the other hand, the 2D case can be thought of as some simplified picture of 
the stagnant layers emerging in liquids that are in contact with a solid body. It is well-known 
(see, e.g. Ref.[71]) that liquids in the close vicinity of a solid interface - at distances within 
a few molecular diameters - do possess completely different physical properties compared to 
those of the bulk phase. In this stagnant region, in which an intrinsically disordered liquid 
phase is spanned by and interacts with the ordering potential of the solid, the liquid's viscosity 
is drastically enhanced and transport processes (related, say, to biased diffusion of charged 
carriers in solutions) are substantially hindered. Thus our model can be viewed as a two-level 
approximate model of this challenging physical system, in which the reservoir mimics the 
bulk fluid phase with very rapid transport, while the adsorbed monolayer (particles on a 2D 
lattice) represents the stagnant layer emerging on the solid-liquid interface. Lastly, the model 
under study represents a novel example [31] of the so-called "dynamic percolation" models 
(see, e.g., Refs. [28-30] and references therein), invoked to describe transport processes in many 
situations with dynamical disorder. In this context, the lattice gas particles can be thought of 
as representing some fluctuating environment, which hinders the motion of an impure molecule 
- the TP, which might be, for example, a charge carrier. A salient feature of this model, which 
makes it different from the previously proposed models of dynamic percolative environments, 
is that we include the hard-core interactions between the "environment" particles and the 
TP, such that the latter may itself influence the dynamics of the environment. As we proceed 
to show, this circumstance turns out to be very important and will result in formation of 
inhomogeneous distribution of the environment particles. 

2.1 Evolution equations. 

Now, we derive the evolution equations in a general, (i-dimensional case. We begin by intro- 
ducing some auxiliary definitions. Let rj = {i](R)} denote the entire set of the occupation 
variables, which defines the instantaneous configuration of the lattice-gas particles at time 
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moment t. Next, let P(R tr ,r/;i) stand for the joint probability of finding at time t the TP 
at the site Rt r and all adsorbed particles in the configuration i]. Then, denoting as rf> v a 
configuration obtained from r\ by the Kawasaki-type exchange of the occupation variables of 
two neighboring sites r and r + e^, and as ff - a configuration obtained from the original 
rj by the replacement ry(r) — > 1 — r/(r), which corresponds to the Glauber- type flip of the 
occupation variable due to the creation/desorption events, we have that the time evolution 
of the configuration probability P(R tr ,ry;t) obeys the following master equation: 

1 d 

d t p(R tr , m t) = — E E { p (Rtr, if";t) - P(Rt r ,m t)} 

H=l r^R tl — e M ,R t r 

+ jE^t 1 " Vfrtr)) P(Rtr - e M , m t)-{l- r/(R tr + e M )) P(R tr , m t)} 
+ ^ E {(l-r/(r))P(Rt r ,f? r ;t)-r ? (r)P(R tr ,r ? ;0} 

+ ^ E {^(r)P(Rtr,^;t)-(l-r / (r))P(R tr ,r ? ;t)}. (8) 

r^Rtr 

The mean velocity V{t) of the TP can be obtained by multiplying both sides of Eq.(8) by 
(Rt r • ei) and summing over all possible configurations (Rt r > f])- This results in the following 
exact equation determining the TP velocity: 

V ^ = J t ^( K tr-e 1 )P(K tr , V ;t) = ^[ Pl (l-k(e 1 ;t^ -p_x(l - fc(e_i;t)) }, (9) 

Rtr,?? 

where 

fc(A; t) = E ^( R tr + A)P(R tr , r?; t) (10) 

Rtr,*? 

is the probability of having at time t a lattice gas particle at position A, defined in the frame of 
reference moving with the TP. In other words, k(X; t) can be thought of as being the particle 
density profile as seen from the moving TP. 

Equation (9) signifies that the velocity of the TP is dependent on the lattice gas particles 
density in the immediate vicinity of the tracer. If the lattice gas is perfectly stirred, or, in 
other words, if k(X;t) = p s everywhere, (which implies immediate decoupling of Rt r and 77), 
one would obtain from Eq.(9) a trivial mean-field result 

V^ = ( Pl -p^)(l- Ps )-, (11) 

r 

which states that the only effect of the medium on the TP dynamics is that its jump time r 
is merely renormalized by a factor (1 — Ps)^ 1 ; (1 — p s )/r defines simply the mean frequency 
of successful jump events. 
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However, the situation appears to be more complicated and, as we proceed to show, k(X; t) 
is different from the equilibrium value p s everywhere, except for |A| — ► oo. This means that 
the TP strongly perturbs the particles distribution on the lattice - it is no longer uniform and 
some non-trivial density profiles emerge. 

Now, in order to calculate the instantaneous mean velocity of the TP we have to determine 
the mean particles density at the neighboring to the TP sites Rt r + e±i, which requires, in 
turn, computation of the density profile k(X; t) for arbitrary A. The latter can be found from 
the master equation (8) by multiplying both sides by r/(Rt r ) and performing the summation 
over all configurations (Rt r ,r?). In doing so, we find that these equations are not closed with 
respect to k(X;t), but are coupled to the third-order correlations, 

T(A, e„; t) = £ r?(R tr + A)r ? (R tr + e M )P(R tr , 77; t) (12) 

In turn, if we proceed further to the third-order correlations, we find that these are coupled 
respectively to the fourth-order correlations. Consequently, in order to compute V(t), one 
faces the problem of solving an infinite hierarchy of coupled equations. 

We resort then to the simplest non-trivial closure of the hierarchy in terms of k(X;t), 
which has been first proposed in Ref.[44], and represent T(A, e„;i) as 



v(Rtr + A)r?(R tr + e M )P(Rt r , m t) 



J2 r?(Rtr + A)P(R tr ,r/;t) r^Rtr + e M )P(R tr , V ; t) 

= k{X;t)k{e„;t), (13) 

Some arguments justifying such an approximation can be found in Ref.[43]. 
Using the approximation in Eq.(13), we obtain 

2dr*d t k(X;t) = Lk(X;t) + 2df, L = ^ A M V M - 4(f + g), (14) 

which holds for all A, except for A = {0, ±ei, e2 . . . , e^}. One the other hand, for these 
special sites A = e u with v = {±1, 2, ... ,d} we find 

2dT*d t k{e„; t) = (L + A u )k(e„; t) + Af, (15) 

where the coefficients are defined by 

A M (t) = l + — p^l-fcCe,.;*)). (16) 
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Note that Eq.(15) represents, from the mathematical point of view, the boundary conditions 
for the general evolution equation (14), imposed on the sites in the immediate vicinity of 
the TP. Equations (14) and (15) together with Eq.(9) thus constitute a closed system of 
equations which suffice computation of all properties of interest, i.e. the time-dependent 
tracer particle velocity and the particle density profiles as seen from the tracer particle. The 
general approach to solution of coupled non-linear Eqs.(9),(14) and (15) has been discussed in 
detail in Ref. [70] . Here we merely note that despite the fact that using the decoupling scheme 
in Eq.(13) we effectively close the system of equations on the level of the pair correlations, 
solution of Eqs.(14) and (15) still poses serious technical difficulties: Namely, these equations 
are non-linear with respect to the TP velocity, which enters the gradient term on the rhs of 
the evolution equations for the pair correlation, and does depend itself on the values of the 
monolayer particles densities in the immediate vicinity of the TP. Solution of this system of 
non-linear equations for several particular case will be presented in the next sections. 



3 Biased tracer diffusion in a ID lattice gas. 

Consider first biased tracer diffusion in the simplest case of a ID hard-core lattice gas with 
conserved particles number; that is, a lattice gas in which all particles were initially introduced 
onto the lattice and particle exchanges with the vapor phase are forbidden (/ = g = 0). Hence, 
the particle density p s on the lattice is kept constant all the time. For simplicity, we also set 
here the spacing a = 1 and r = r* = 1. These parameters can be trivially restored in final 
results. 

In absence of external field acting on the TP, when its transition probabilities are sym- 
metric, tracer dynamics has been studied first in Refs. [18-22] and the TP mean-square dis- 
placement < X 2 (t) > has been determined exactly: 

< xHt) >_ (IzfJ M (17) 

Ps V 7T 

which thus follows an anomalous diffusion law in Eq.(2) with a = 1/2. Consequently, the 
" diffusivity" D(t) is given here by 

D(t) = ( 1 ~ ps ^ * ^ as i^oo. (18) 
Ps V2rt 

Now, if we charge the TP and switch on an external electric field, how will the TP 
velocity and the mobility behave? An early incomplete answer was presented in Ref. [72] 
which focused on the extreme case of an infinitely strong electric field E = +oo such that 



10 



the TP performs a totally directed random walk constrained by hard-core particles. By 
noticing that such a model is tantamount to the ID Stefan freezing problem or the problem 
of directional solidification [73], (the role of the latent heat being played by the lattice gas 
particles), it was found that in this extreme case the TP mean velocity follows 

V(t) = ^, (19) 

where 7^ is a function of p s defined implicitly by 

Vexp(^)erfc(^|) = l-p s , (20) 

erfc(x) being the complementary error function. Hence, the TP mean velocity in the biased 
case decays at exactly the same rate as the TP " diffusivity" in the unbiased case! In view of 
such a coincidence, one is, of course, prompted to consider the general case E < 00 and to 
check the prefactors in Eq.(19) having in mind to verify the validity of Eq.(5) for anomalous 
tracer diffusion in such random dynamical environments [74]. 

The general case of arbitrary < E < 00 has been studied in Ref . [44] within the framework 
of the approach outlined in Section 2 and subsequently, in terms of a rigorous probabilistic 
approach in Ref. [46]. Results of these two approaches coincide, which justifies the decoupling 
approximation in Eq.(13). It has been found that in this general case the TP mean velocity 
V(t) obeys at sufficiently large times the form of Eq.(19), i.e. V(t) = 7_e/v / ^ in which law 
the prefactor 7^ is determined by the following equation [44,46]: 

i+M - 1 + p 



l-exp(-0E) i 
T u, a 1 1 1 Ps \ __ exp(/3£)p2 



where 



In the limit (3E <C 1 the solution of the transcedental Eq.(21) has the form 



(22) 



^ (l-Ps) [2(3E 
1E ~^^^ — ' (23) 



which implies that the time-dependent mobility in this case is given by 



m = S^4=. (24) 

Ps V2lTt 
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Finally, by comparing the prefactors in the latter equation and Eq.(18), we find that the 
generalized Einstein relation in Eq.(5) holds exactly for such an anomalous diffusion! 

To close this section, we make several comments on the particle density distribution as seen 
from the tracer particle. Since the TP has a non-zero mean displacement, while the lattice 
gas particles have symmetric hopping probabilities and hence, a zero mean displacement, they 
tend to accumulate in front of the TP. Consequently, a condensed, traffic jam- like region is 
formed in front of the TP, and a depleted by the lattice-gas particles region past the TP 
emerges. As a matter of fact, such a inhomogeneous density distribution is not stationary: 
the size of the traffic jam-like region, in which the particle density profile k(X, t) is nearly 
constant and larger than p s , k(X,t) = p s /(l + I+^e)), grows in proportion to the TP mean 
displacement. The size of the depleted region with a constant density k(X,t) = + 
I-(1e)) < Ps also grows in proportion to < X(t) >= V(t)t. 

4 "Shock" propagation in a non-equilibrium ID lattice gas. 

In this section we address the question of the validity of the generalized Einstein relation in 
a somewhat exotic situation in which the initial particle distributions from both sides of the 
TP has an S'-shape, shock-like profile characterized by two different mean densities [45]. That 
is, as depicted in Fig.2, we suppose that the TP is initially at the origin, and the particles 
mean densities from the left and from the right of the origin, which we denote as p_ and p + , 
respectively, are not equal to each other. The TP thus defines position of the "shock front". 
Without of lack of generality, we suppose that P- > p+ > 0, and will call in what follows 
the particle phase which initially occupied the left half- line as the high-density phase (HDP), 
while the phase initially occupying the right half-line will be referred to as the low-density 
phase (LDP). 

Further on, after deposition on the lattice, particles start to perform symmetric, (i.e. 
with equal probabilities for going to the left or to the right), hopping motion between the 
nearest lattice sites under the constraint that neither two particles can simultaneously occupy 
the same lattice site and can not pass through each other. Next, we stipulate that the TP 
is subject to a constant force F which favors its jumps in a preferential direction. Thus for 
the TP the probabilities of going to the right (p) and to the left (q) will be different from 
each other. We will suppose that F is directed towards the HDP. We adopt the convention 
that in this case F is negative, F < 0; the TP hopping probabilities p and q are related to 
the force and the reciprocal temperature f3 through p/q = exp((3F) and p + q = 1. From 
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The "shock front" particle 
p_ I Pt- 

1/2 1/2 q P 1/2 1/2 
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Figure 2: Initial " shock" configuration of the lattice gas particles on a ID lattice. p_ is the initial 
mean density of particles at the half-line — oo < X < and /?+ is the mean density of particles at the 
half-line < X < oo. All lattice gas particles, excluding the TP (the "shock" front particle), have 
symmetric transition probabilities. The TP has asymmetric transition probabilities and experiences 
the bias in the negative direction. 

the physical point of view, such a constant force can be understood as an effective boundary 
tension derived from the solid-on-solid-model Hamiltonian of the phase-separating boundary 
[75-77] and mimic, in a mean-field fashion, the presence of attractive interactions between 
the lattice-gas particles which are not explicitly included in the model. Relation between p, q 
and the amplitude of attractive interparticle interactions has been discussed in Refs.[78] and 
[79]. 

The TP dynamics in such a system has been analysed in Refs.[45] and [46] and it has been 
found that the TP mean velocity obeys < V(t) >= ^p/yt at sufficiently large times, where 
the prefactor "fp is determined implicitly as the solution of the following equation: 

(IP- V P+ fn <,\ 

- q - p, (25) 



1 + I_(7f) 1-/+(7f) 
which generalizes the result of the previous section for the case of p_ 7^ p+. 

Equation (25) predicts that three different regimes can take place depending on the relation 
between p/q and p±: 

(1) When p(l — p + ) > q(l — p_) the parameter is finite and positive, which means 
that the HDP expands compressing the LDP. In the particular case p/q = 1 and p + = the 
parameter 7^ appears to be a positive, logarithmically growing with time function, which 
behavior agrees with the results of Ref.[80]. 

(2) When p(l — p + ) < q(l — p_) the parameter is less than zero - the expanding LDP 
and the applied force effectively compress the HDP. 

(3) When p(l — p + ) = q(l — p_), the parameter = 0. This relation between the system 
parameters when the HDP and the LDP are in equilibrium with each other and the TP mean 
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displacement is zero, was found also in Ref.[42] from the analysis of the stationary behavior 
in a finite ID lattice gas. 

"Shock front", or, in other words, the TP propagation in cases (1) and (2) has been 
discussed in detail in Ref.[44]. Here, in view of our interest in the Einstein relation, we will 
focus on the case (3) in which situation the high- and low-density phases coexist such that 
the TP mean displacement is equal to zero. On the other hand, it is clear that also in this 
case the TP will wander randomly around the equilibrium position and one may expect that 
its mean-square displacement will not be equal to zero. 

To determine the TP mean-square displacement < X 2 (t) > we may try to extend the 
approach outlined in Section 2 which will require some more cumbersome analysis. We will 
however pursue a different, heuristic approach and instead of evaluating < X 2 {t) > in terms 
of the approach of Section 2 we will resort to the Einstein relation in Eq.(5) assuming that 
it also holds in the inhomogeneous situation under study. Our result will be then checked by 
numeric simulations. 

To do this, we note first that contrary to the model of Section 3, here the TP does not 
move, on average, when the external force F acting on this particle is not equal to zero, but 
rather to some critical value F c , such that 

^'Kw;)- < 26 > 

Hence, it seems reasonable to define the TP mobility in this inhomogeneous case as 

fi(t) = lim 77^%- = 4= lim TTT^TTT (27) 
m ' f->f c (F - F c ) y/t F-^Fc (F -F c ) y ' 

Further on, we find from Eq.(25) that in the limit 0(F — F c ) <C 1 the prefactor 7^ follows 

/ 2 p(l-f+)-?(l 
vr qp- + pp + 



which implies that p(t) obeys 



„<«)-/» Z 1 -"->('-"+> JI. (29 , 

Now, assuming that Eq.(5) holds, we get 

<*(«)>= . (l-'-H 1 -'^ JE, (30, 
(P- + P+ ~ 2 p- p+) V vr 

which reduces to the classical result in Eq.(17) in the limit p_ = p + . 

In Fig. 3 we present results of extensive Monte Carlo simulations which confirm the result 
in Eq.(30) and hence, signify that the Einstein relation seemingly works in this inhomogeneous 
case. 



14 












20 



40 



1/2 



60 



80 



100 



t 



Figure 3: Mean-square displacement of the TP in the critical case p/q = (1 — p~)/(l — p+)- The 
solid lines show our analytical prediction from Eq.(30) and the noisy lines present the results of 
Monte Carlo simulations. The curves from top to bottom correspond to the following values of the 
parameters (p_,p + ): The first two curves are the analytical result in Eq.(17) and numerical data 
for the symmetric cases (0.4,0.4) and (0.5,0.5). The lower curves correspond to (0.6,0.5), (0.9,0.4), 
(0.7,0.5) and (0.8,0.2) respectively. 

5 Biased TP diffusion in a ID lattice gas in contact with a 
reservoir. 

We turn now to discussion of the TP biased diffusion in situation in which the reservoir 
is present and particle exchanges with the reservoir are allowed; that is, when both the 
parameters / and g are not equal to zero. 

In ID systems with continuous particles exchanges with a reservoir, which situation is 
appropriate to adsorption on polymer chains [81,82], the biased TP will ultimately move 
with a constant velocity V and particles distribution around the TP will be characterized by 
stationary density profiles k(X) [43]. Here, the general solution of Eqs.(14) and (15) has the 
following form: 




(31) 



where the characteristic lengths A± obey 



A± = =F o l n 



A 1 + A_i + 2(/ + g) T J (A, + A_i + 2( f + gj) 



-AAiA-x 



(32) 



2A 1 
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while the amplitudes K± are given respectively by 

A 1 - A-i A x - A_ x 

K +=Ps- A -, ; rr^, and K- = p s - — — — . (33) 

v4_i - Ai exp(-cr/A + ) A-i exp(-cr/A_) - A 1 

Note that A_ > A + , and consequently, the local density past the TP approaches its non- 
perturbed value p s slower than in front of it; this signifies that correlations between the TP 
position and particle distribution are stronger past the TP. Next, K + is always positive, 
while K_ < 0; this means that the density profile is a non-monotoneous function of A and is 
characterized by a jammed region in front of the TP, in which the local density is higher than 
p s , and a depleted region past the TP in which the density is lower than p s . This is, of course, 
quite similar to the behavior observed in situations in which exchanges with the reservoir are 
not allowed (Sections 3 and 4); essential difference is, however, that here the inhomogeneous 
particles distribution around the TP attains a stationary form. 

0.20 
V 

0.15 
0.10 
0.05 
0.00 

0.0 0.5 1.0 

/ 

Figure 4: Terminal velocity of the TP as a function of the "creation" probability f at different values 
of the parameter g. The TP hopping probabilities are p+i — 0.6 and p_i = 0.4. The solid lines give 
the analytical solution while the filled squares denote the results of Monte-Carlo simulations. Upper 
curves correspond to g — 0.8, the intermediate - to g — 0.5 and the lower - to g = 0.3, respectively. 

Now, we are in position to obtain a system of two closed-form non-linear equations deter- 
mining implicitly the unknown parameters Ai and A-i, which will allow us to compute the 
TP terminal velocity, related to A±i through V = a(A 1 - A_i)/2t* . Substituting Eq.(31) 
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into Eq.(16), we find 



1 - p s - p s 



and 



A_iexp(o-/A' + ) - Ai 

I-Ps-P. 



(34) 



(35) 



A_i - A 1 exp(a/X'_) 

which system of equations defines the stationary velocity V of the TP (see Fig.4) as well as 
the density profiles (see Fig.5). For arbitrary values of p, f and g the parameters A±i, defined 
by Eqs.(34) and (35), and consequently, V can be determined only numerically (see Figs. 4 to 
6). However, V can be found analytically in the explicit form in the limit of a vanishingly 
small force E, [3E <C 1. Expanding A±i in the Taylor series in powers of E and retaining 
only linear with E terms, we find that the TP velocity follows 

V ~ C'E, (36) 

which relation can be thought off as the analog of the Stokes formula for driven motion in 
a ID lattice gas undergoing continuous particles exchanges with the vapor phase. Eq.(36) 
signifies that the frictional force exerted on the TP by the lattice gas particles is viscous. 
The friction coefficient, i.e. the proportionality factor in Eq.(36), is given explicitly by 



2r 



C = 11 i + PsT " (37) 

Pa\l-p s )V T (f + g) 1 + ^1 + 2(1 + r*(l - Ps)/T)/{f + g) J 

Note that ( in Eq.(37) can be written down as the sum, ( = ( cm + ( COO p- The first term, 
( cm = 2r//3cr 2 (l— p s ) is a typical mean-field result and corresponds to a perfertly homogeneous 
lattice gas (see discussion following Eq.(lO)). The second term, 

8r*p s 1 
Cc ° OP = (3^(1 - Ps )(f + g)i + yiT2(l + r*(l - p.)/T)KF+gl' { ' 

has, however, a more complicated origin. Namely, it reflects a cooperative behavior emerging 
in the lattice gas, associated with the formation of inhomogeneous density profiles (see Fig.5) 
- the formation of a "traffic jam" in front of the TP and a "depleted" region past the TP 
(for more details, see Ref.[43]). The characteristic lengths of these two regions as well as 
the amplitudes K± depend on the magnitude of the TP velocity; on the other hand, the 
TP velocity is itself dependent on the density profiles, in virtue of Eq.(9). This results in 
an intricate interplay between the jamming effect of the TP and smoothening of the created 
inhomogeneities by diffusive processes. Note also that cooperative behavior becomes most 
prominent in the conserved particle number limit [44]. Setting /, g — > 0, while keeping their 
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ratio fixed (which insures that p s stays constant), one notices that Ccoop gets infinitely large. 
As a matter of fact, as we have shown already in section 3, in such a situation no stationary 
density profiles around the TP exist. 

0.5 r 

k 

0.4 - 

0.3 - 

0.2 - 

0.1 - 

-15 -10 -5 5 10 15 

n 

Figure 5: Density profile around stationary moving TP for f — 0.1, g — 0.3 and p = 0.98. The solid 
line is the plot of the analytical solution. Filled squares are the results of Monte-Carlo simulations. 

In order to check our analytical predictions, we have performed numerical Monte Carlo 
simulations. Results of these simulations, performed at different values of the parameters /, 
g, and pi, are also represented in Figs. 4 and 5. 

Consider finally the situation with E = 0, in which case the terminal velocity vanishes 
and one expects conventional diffusive motion with the mean square displacement of the form 
of Eq.(l). Heuristically, we can evaluate D for the system under study if we assume the 
validity of the Einstein relation in Eq.(5). Noticing that here the biased TP mobility p is just 
p = 1/C, we find [43]: 

D _ a 2 {l- Ps ) f i Ps t* 2 

~ 27 \ + r(f + g) 1 + ^1 + 2(1 + r*(l - p s )/r)/(f + g) 

Monte Carlo simulations (see Fig. 6) evidently confirm our prediction for D given by Eq.(39), 
and hence, confirm the validity of the Einstein relation for the system under study. This is, 
of course, not an unexpected, but still a non-trivial result [42]. 
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Figure 6: Self-diffusion coefficient of the TP as a function of the adsorption probability f. Notations 
and values of g are the same as in Figs. 4. 

6 Biased TP diffusion in 2D and 3D lattice gases. 

We turn now to the analysis of a biased TP diffusion in higher dimensional (2D and 3D) lattice 
gases exposed to a vapor phase. Here, the situation gets somewhat more difficult from the 
computational point of view; we have now to solve the partial difference equations problem 
(14), (15) rather than the mere difference equations arising in the ID case. 

Below we briefly outline such a solution for a 2D case. Results for three-dimensions 
are obtained along essentially the same lines and here we display them without derivation. 
Solution in two-dimensions can be found in a most convenient fashion if we introduce the 
generating function for the particle density profiles, defined as 



Multiplying both sides of Eqs. (14) and (15) by w^w^ 2 , and performing summations over n\ 
and ri2, we find that H(wi,W2) is given explicitly by 




(40) 



«i=- oo ri2 =— oo 




(41) 



where a = Y^ u A v + 4(/ + g) and 



K(w u w 2 ) = Y,Mw^ Ul ~ 1) (*(e„) - Ps) + Ps{A l - A_i)(u;i - w^ 1 ). 



(42) 
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Equations (41) and (42) determine the generation function for the density profiles exactly. 

Before we proceed to the inversion of H(w\,W2) with respect to the variables w\ and u>2, 
we note that we can access interesting integral characteristics of the density profiles directly 
using the result in Eqs.(41) and (42). Namely, as we have already remarked, the presence of 
the driven TP induces an inhomogeneous density distribution in the lattice. One can thus 
pose a natural question whether equilibrium between adsorption and desorption processes 
gets shifted due to such a perturbancy, i.e. whether the equilibrium density in the lattice 
gas is different from that given by Eq.(6). The answer is trivially "no" in the case when the 
particles number is explicitly conserved, but in the general case with arbitrary / and g this 
is not at all evident: similarly to the behavior in ID system one expects that also in two- 
dimensions the density profiles are asymmetric as seen from the stationary moving TP and 
are characterized by a condensed, "traffic-jam" -like region in front of and a depleted region 
past the TP. One anticipates then that the desorption events are favored in front of the TP, 
while the adsorption events are evidently suppressed by the excess density. On the other 
hand, past the TP desorption is diminished due to the particles depletion while adsorption 
may proceed more readily due to the same reason. It is thus not at all clear a priori whether 
these two effects can compensate each other exactly, in view of a possible asymmetry of the 
density profiles, as it happens in the ID model (see Fig.5). 

For this purpose, we study the behavior of the integral deviation of the density from the 
equilibrium value p s , i.e. 17 = ^ni°^-<x> Sn 2 °^-oo (p* ~ M n i) n 2)), which can be computed 
straightforwardly from Eqs.(41) and (42) by setting both w\ and W2 equal to unity. Noticing 
that K(wi = 1, W2 = 1) = 0, and that A\ + A-i + 2A2 — a = — 4(f + g), we obtain then that 
Jl is strictly equal to 0. This implies, in turn, that the perturbancy of the density distribution 
in the lattice gas created by the driven TP does not shift the global balance between the 
adsorption and desorption events. 

Inversion of the generating function with respect to w\ and u>2 requires quite an involved 
mathematical analysis, which has been presented in detail in Ref.[70]. General solution for 
the density profiles reads: 




Ps 



+ 




V 




(43) 



with 



) 



ni/2 



/ 

Jo 




F, 



n\,ni — 



A 2 t) dt 



(44) 
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where l n (z) stands for the modified Bessel function. We mention that F nim has an interesting 
physical interpretation in terms of the generating function of a random walk of a single particle 
(that is, in absence of the particles environment) [83]. 

Now, the Eqs.(43) and (44) display k(ni,n 2 ) as a function of the coefficients A v that 
remain to be determined. As a matter of fact, these coefficients depend themselves on the 
local densities in the immediate vicinity of the tracer, i.e. on k{e v ). This implies that we 
have to determine them from Eqs.(43) and (44) in a self-consistent way [70]. Some analysis 
(see Ref . [70] ) shows that A v are determined implicitly as the solution of the following system 
of three non-linear matrix equations 



Vi/ = {±1,2}, A h 




p s - p s (A 1 - 



detCv 
detC 



(45) 



where 



C = 



(A, 



AiV-iFe^ 



a 



A 2 V^ 2 F ei \ 
A 2 V- 2 F e _ 1 



(46) 



\ ^iV_iF e2 A_!V!F e2 A 2 V„ 2 F e2 -a/ 

the matrix C„ stands for the matrix obtained from C by replacing the v-th column by the 
column- vector F, 

I (Vi - V_i)F ei \ 
F= (Vi-V-i)^ , (47) 
V (Vx - V_i)F e2 / 

while fe(e„) are expressed in terms of A v as 



1 + 



Ar*p u 



Lastly, we find that in 2D the TP terminal velocity obeys 

4t* pt det Ci - p-! det C_i i - 1 



V = -(pi-p-l)(l-ps){l+Ps — 

t L r 



(48) 



(49) 



detC J 

which represents the general force-velocity relation for the system under study, valid for 
arbitrary magnitude of the external bias and arbitrary values of other system's parameters, 
except apparently the situation when r* = oo, /, g = and their ratio f/g = const. This 
case corresponds to the percolation-type problem describing biased diffusion in a randomly 
diluted lattice - the model which has been discussed in the Introduction. Some analysis shows 
that V defined by Eq.(49) does not vanish for any value of external bias E. This signifies that 
the decoupling approximation underlying the derivation of the result in Eq.(49) is, of course, 
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inappropriate for situations with quenched disorder, but only for situations in which some 
mixing, either due to exchanges with the reservoir, or due to diffusional processes, is present. 

We turn now to the limit [3E <C 1, in which case the problem simplifies considerably and 
allows to obtain explicit results for the local densities in the immediate vicinity of the TP 
and consequently, for the TP terminal velocity and its diffusivity. 

In this limit, we arrive again at a Stokes-type formula of the form V ~ E/(, where now 



<- 4T 



with 

1 l + r*(l-p a )/r 
x = -- 



l + T - 



(/ + 5 + l + r*(l-p a )/r) [c 2 {x)-x) ^ 



(50) 



2f + g + l + T*{l- Ps )/t' 



and C 2 {x) = iyJ^ e~* ((Io(zi) - I 2 (xt))l (xt)dt 



(51) 

Note that we again are able to single out two physically meaningful contributions to the 
friction coefficient £• Namely, the first term on the rhs of Eq.(50) is just the mean-field- type 
result corresponding to a perfectly stirred lattice gas, in which correlations between the TP 
and the lattice gas particles are discarded. The second term mirrors the cooperative behavior 
emerging in the lattice gas and is associated with the backflow effects. In contrast to the 
ID case, however, the contribution to the overall friction coefficient stemming out of the 
cooperative effects remains finite in the conserved particles number limit. 

Lastly, assuming a priori that the Einstein relation holds for the system under study, we 
find the TP diffusion coefficient: 

D = " Ps) |l " r 4{f + g + 1 + T * (1 _ Ps )/ T )C 2 (x) - 1 + (3 Ps - l)r*/r J ' (52) 
It seems now interesting to compare our general result in Eq.(52) against the classical result 
of Nakazato and Kitahara [60], which describes TP diffusion coefficient in a 2D lattice gas 
with conserved particles number. Setting / and g equal to zero, while assuming that their 
ratio has a fixed value, f /g = p s j (1 — p s ), we find that the right-hand-side of Eq.(52) attains 
the form 

4r U Ps) \ L t l + (l-p s )rVr-(l-2/vr)(l + (l-3p s )rVr)/' ^ 
which expression coincides with the earlier result obtained in Refs.[60] and [64] within the 
framework of a different, compared to ours, analytical technique. The result in Eq.(53) is 
known to be exact in the limits p s <C 1 and p s ~ 1, and serves as a very good approximation 
for the self-diffusion coefficient in hard-core lattice gases of arbitrary density [61]. 
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In a similar fashion, we find that in 3D the TP terminal velocity obeys V ~ E/(, where 
/3a (I - p s ) p a 2 l{l _ p s ).( ao c 3 (2A /a ) - A j 

with 

A) = 1 + ^(1 - Ps), ao = 6(lH — 1 — J, (55) 

and 

C 3 (x) = y^ e-*l2(xt)(I (xt)-I 2 (xt))dt| ={P(0;3x)-P(2e 1 ;3x)}- 1 , (56) 
P(R;£) being the generating function, 

+oo 

P(R;0 = J2 P ^^ ( 57 ) 

3=0 

of the probability Pj (R) that a walker starting at the origin and performing a Polya random 
walk on the sites of a 3D cubic lattice will arrive on the j-th step to the site with the lattice 
vector R [83]. 

Assuming next that the Einstein relation holds, we find that in 3D the TP diffusion 
coefficient reads: 

D = <^> j, _ ^f {aoCl(2Aa/aa) _ , + 1^2)-'} ,58) 

In the particular case of conserved particles number, when /, g — > but their ratio f/g is 
kept fixed, f/g = — p s ), the latter equation reduces to 



(59) 



which coincides again with the classic result obtained earlier in Refs.[60] and [64]. 

Asymptotical behavior of the density profiles at large distances from the TP follows from 
the analysis of the analyticity properties of the complex function N(z) = Y2n=-oo k(n, 0)z n . 
It has been shown in Ref.[70] that in front of the TP, k(n, 0) always decays exponentially with 
the distance: 

exp ( — n/A H 



k(n,0)~p s + K + ^— (60) 

where the characteristic length A + obeys, in two-dimensions: 

^ = ^~ 1 (t-\^- a 2 + \H^-A2)'-a 1 a.A) 
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Note that A + stays finite for any values of the system parameters. On contrary, the behavior 
of the density profiles at large distances past the tracer qualitatively depends on the physical 
situation studied. In the general case when exchanges with the particles reservoir are allowed, 
the decay of the density profiles is still exponential: 

exp ( — n/A_ ) 

fc(-n,0)~p s + K_ ^ J -, (62) 

where 

*-=-*-' (£{f-W(!- j4 ! '- Aij1 - 1 }) (63) 

Note that in the general case the characteristic lenghts again, similarly to the ID case, satisfy 
the inequality A_ > A+, which means that the correlations between the TP and the lattice 
particles are always stronger past the TP than in front of it. 

Such correlations become extremely strong in the special case when the particles exchanges 
with the vapor phase are forbidden, i.e., in the conserved particles number limit. In this case, 
we have that A_ becomes infinitely large and, in the limit n — > +00, the particle density 
follows 

fc(-n,0)=p s -^^(l + |- + o(-L)), K'_(d = 2)>0. (64) 

Remarkably enough, in this case the correlations between the TP position and the particles 
distribution vanish algebraically slow with the distance! This implies, in turn, that in the 
conserved particles number case, the diffusive mixing of the lattice gas is not efficient enough 
to prevent the appearence of the quasi-long-range order and the medium "remembers" the 
passage of the TP on a long time and space scale, which signifies very strong memory effects. 
In three dimensions, an analogous result for the asymptotical behavior of the density profiles 
as seen from the stationary moving TP is: 

K'_ (d = 3) ln(n) , s 
fc(-n,0,0)~p s -i ' V \ K^(d = 3)>0, (65) 

where K'_(d = 3) is a constant [31]. 

7 Biased, single vacancy-mediated TP diffusion. 

We turn finally to the extreme case of a very dense lattice gas in which all sites except one 
(called in what follows "a vacancy") are filled with identical neutral, hard-core particles (see, 
Fig. 7). The particles move randomly by exchanging their positions with this single vacancy, 
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subject to the hard-core exclusion constraint. The just described model, which represents, in 
fact, one of the simplest cases of the so-called "slaved diffusion processes", has been studied 
over the years in various guises, such as, for instance, the "constrained dynamics" model of 
Palmer [85]. 

E 

► 




Figure 7: Two-dimensional, infinite in both directions, square lattice in which all sites except one are 
filled with identical hard-core particles (gray spheres ) . The black sphere denotes a single tracer particle, 
which is subject to external held E, oriented in the positive x\ direction, and thus has asymmetric 
hopping probabilities. The arrows of different size depict schematically the hopping probabilities; a 
larger arrow near the TP indicates that it has a preference for moving in the direction of the applied 
held. 

Brummelhuis and Hilhorst [62] were first to present an exact solution of this model in 
the lattice formulation. It has been shown that in the presence of a single vacancy the TP 
trajectories are remarkably confined; the mean-square displacement shows an unbounded 
growth, but it grows only logarithmically with time, 

2 ln(n) 

< R n >~ — ; — , as n — ► oo, (66) 

TTyK — 1) 

which implies that the TP diffusivity D n obeys 

ln(n) 
47r(7r — l)n 

Brummelhuis and Hilhorst have also found [62] that at sufficiently large times the probability 
P n (R n ) of finding the TP at time moment n at position R n converges to a limiting form as 



a function of the scaling variable r\ = |R n |/y / ln(n). Still striking, this limiting distribution is 
not a Gaussian but a modified Bessel function Ko(rj), which signifies that the successive steps 
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of the TP, although separated by long time intervals, are effectively correlated. These results 
have been subsequently reproduced by means of different analytical techniques in Refs.[86, 87] 
and [88]. 

As in the previous sections, we focus here on TP dynamics in the biased case, when the TP 
is charged, while the lattice gas particles are neutral, and the system is subject to a constant 
electric field. We aim to determine exactly the TP mean-displacement < R„ >, which will 
give us an access to the TP mobility and thus will allow to verify the validity of the Einstein 
relation in this model with remarkably confined dynamics. 

A standard approach to define the properties of the TP random walk, as exposed in 
Section 2, would be to start with a master equation determining the evolution of the whole 
configuration of particles. In doing so, similarly to the analysis of the tracer diffusion on 
2D lattices in the presence of a finite vacancy concentration (see, e.g., Ref.[43]), one obtains 
the evolution of the joint distribution P n (R n , Y) of the TP position R„ and of the vacancy 
position Y at time moment n. The property of interest, i.e. the reduced distribution function 
of the TP alone will then be found from P n (R n , Y) by performing lattice summation over all 
possible values of the variable Y. 

Here we pursue, however, a different approach, which has been first put forward in the 
original work of Brummelhuis and Hilhorst [62]; that is, we construct the distribution function 
of the TP position at time n directly in terms of the return probabilities of the random walk 
performed by the vacancy. The only complication, compared to the unbiased case considered 
by Brummelhuis and Hilhorst [62], is that in our case ten different return probabilities will be 
involved, instead of three different ones appearing in the unbiased case. Hence, the analysis 
will be slightly more involved. 

We begin by introducing some basic notations. Let 

• P ra (R n ) be the probability that the TP, which starts its random walk at the origin, 
appears at the site R n at time moment n, given that the vacancy is initially at site Yo. 

• F*(0 | Yo) be the probability that the vacancy, which starts its random walk at the 
site Yo, arrives at the origin for the first time at the time step n. 

• F*(0 | e„ | Yo) be the conditional probability that the vacancy, which starts its random 
walk at the site Yo, appears at the origin for the first time at the time step n, being at 
time moment n — 1 at the site e v . 
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Further on, for any time-dependent quantity L n we define the generating function: 

+00 



L(0 = J2 L nC, (68) 

n=0 

and for any space-dependent quantity Z(X) the discrete Fourier transform: 

Z(k) = £exp(i(k-X))y(X), (69) 
x 

where the sum with the subscript X = (xi,x 2 ) runs over all lattice sites. 

Now, following Brummelhuis and Hilhorst [62], we write down directly the equation obeyed 
by the reduced probability distribution P n (R n ) (cf Ref.[87] for a study of the joint probability 
of the TP position and of the vacancy position in the unbiased case): 



Pn(Rn) = <5 R „,0 \^-J2F*(0\Y )J + 

+00 +00 +00 +00 

+ ••• < W+-+m p+ i,n ^ ■ ■ ■ yi ^1 +-+e„ p ,Rn x 

p=lmi=l m p =lm p +i=0 v\ v v 

m p+l 

1- £i7(0|-e„ p ) 

j=0 

x ^(0|e„ p | -e I/p _ 1 )...F^(0|e I/a | - e„ x )F* mi (0\e^ | Y„). (70) 

Next, using the definition of the generating functions and of the discrete Fourier transforms, 
Eqs.(68) and (69), we obtain the following matricial representation of the generating function 
of the TP probability distribution: 



P(k;0 = ^ ^l + p-\k;0^^(k;0^*(0 | e M | Y ;^ • 
In Eq.(71) the function P(k; £) stands for the determinant of the following 4x4 matrix, 



(71) 



P(k;0 = det(l-T(k;£)), 
where T(k; £) is a 4 x 4 matrix of the form: 



(72) 



T(k;£) = 



( e ifcl Ai_i(0 e*A 1;1 (£) e ifcl Ai,_ 2 (0 e*A li2 (£) \ 

r ifc iA_i_i(0 e - ifcl A_i,i(£) e-*A_ 1 ,_ 2 (e) e- ifcl A_i )2 (0 
e ita A 2 _i(0 e ifc2 A 2>1 (0 e ifc2 A 2 ,_ 2 (0 e lfe2 A 2 , 2 (0 

Ve- ifc »A_ 2 _i(0 e-^A-a,!^) e-^A_ 2 ,_ 2 (0 e~^A^ 2 (0j 



(73) 
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and the coefficients A,^(£), v,n = ±1,±2, defined as 

+00 

A^(0 = F*(0 I e„ I e M ;0 = J>n*(0 | e„ \ e M )f\ (74) 

n=0 

are the generating functions of the conditional probabilities for the first time visit of the origin 
by the vacancy. Lastly, the matrix f/ M (k;£) in Eq.(71) is given by 

U,(k; = V(k; £(1 - e-( k ^))(I - T(k; 0)u,l e^-l (75) 

V 

The property of interest - the TP probability distribution function, will be then obtained by 
inverting P(k; £) with respect to the wave-vector k and to the variable £: 

p »( R »» - hi i^w?- /> £ dfe ^'^o- < 76 » 

where the contour of integration C encircles the origin counterclockwise. 

After straightforward but rather tedious calculations (see, Ref.[90]), we find that in the 
small-k limit and £ — ► 1~, which defines the leading large-n behavior of the probability 
distribution, the generating function P(k;£) obeys 



P(k;0 = T^j 1 " (-*aofci + ^aifc? + ^a 2 fci)ln(l-0} \ 



(77) 



where the coefficients 

'<*o(£) = 7r- 1 sinh(/3P/2)((2^ - 3)cosh(/3P/2) + l) -1 , 
ai(P) = vr- 1 cosh(/3P/2)((2vr - 3) cosh(/3P/2) + l)" 1 , (78) 
a 2 (P) = vr' 1 (cosh(/3P/2) + 2vr - 3) _1 , 

are all functions of the field strength E and of the temperature only. 

Now, the leading large-n asymptotical behavior of the TP mean displacement can be 
obtained directly from Eq.(77), since the generating function of the TP mean displacement: 

+00 

< R(0 >^ 2 < R « > ( 79 ) 



n=0 



obeys (see, e.g., Ref.[61]): 



< R «)>=-' £ (0;«e 1+ |£(0; f )e 2 . (80) 
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Consequently, differentiating the expression on the right-hand-side side of Eq.(77) with re- 
spect to the components of the wave-vector k, we find that the asymptotical behavior of the 
generating function of the TP mean displacement in the vicinity of £ = 1~ follows 

Further on, using the discrete Tauberian theorem (cf. Ref.[83]) and Eq.(78), we find the 
following force- velocity relation for the system under study 

^ ( / i\ i \ /I sinh(/?£/2) , \ 

< R„ >~ (a (E) Inn) e, = (- - _ - + , Inn) e 1; as n ^ oo, (82) 

which shows that the TP mean displacement grows logarithmically with n. In the limit 
f3E <C 1, the coefficient ao(E) obeys 

a (E) = 1 -^ T)+ O(E% (83) 

and hence, the mobility /j, n , defined in Eq.(3), follows 

P In (n) 

V™ ~ Tl TT ' as n ^ oo. (84) 

Att[jt — 1) n 

Comparing next the result in Eq.(84) with that for the diffusivity D n , Eq.(67), derived by 
Brummelhuis and Hilhorst [62] for the unbiased case, we infer that the TP mobility and 
diffusivity do obey, at least in the leading in n order, the generalized Einstein relation in 
Eq.(5)! Note, that this can not be, of course, an a priori expected result, in view of an 
intricate nature of the random walks involved and anomalous, logarithmic confinement of 
the random walk trajectories. 

We turn next to calculation of the asymptotic forms of the probability distribution P n (R n ). 
Inverting P(k;£) with respect to k, we find that in the large-n and large-X limits, in the 
general case when > E < oo, P n (R n ) obeys 

P n (R n ) ~ (Try/ ai (E)a 2 {E)\n (n)) exp (^|f|^i) K ( VE (n)), (85) 
where Kq is the modified Bessel (McDonald) function of zeroth order, while 



/ 2 c£(E) / x? xl 

^ s i^ + ^^W) + ^(ify (86) 

x\ and X2 being the components of the vector R n . The distribution in Eq.(85) reduces to the 
result obtained by Brummelhuis and Hilhorst [62] in the limit E = 0. 
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8 Conclusions. 



To conclude, in this presentation we have overviewed some recent results on biased tracer 
diffusion in disordered dynamical environments represented by hard-core lattice gases. We 
have considered several possible situations including ID lattice gases with conserved and non- 
conserved particles density with initial homogeneous or "shock" -like particle distributions, 
d-dimensional lattice gases undergoing continuous exchanges with a reservoir, as well as an 
extreme case of biased tracer diffusion in a 2D highly-packed lattice gas containing only a 
single vacancy. The main emphasis of our analysis here has been put on the question whether 
the Einstein relation between the tracer particle mobility and the diffusivity is valid, despite 
the fact that in some cases the tracer diffusion is anomalously confined or, in case of regular 
diffusion, the mobility and the diffusivity are non-trivial functions of the lattice gas density 
and other pertinent parameters. For some situations we were able to furnish exact solutions, 
which show explicitly that this fundamental relation holds. In other instances, to calculate 
the tracer particle mobility we invoked an approximate approach based on decoupling of 
the third-order correlation functions. The tracer particle diffusivity has been then obtained 
by assuming that the Einstein relations holds. Extensive numerical simulations have been 
performed, which confirm our analytical predictions and thus the validity of the Einstein 
relation. 

Apart of it, we have shown that the lattice gas particles distribution as seen from moving 
TP is strongly inhomogeneous: the local particle density in front of the TP is higher than the 
average, which means that the lattice gas particles tend to accumulate in front of the driven 
TP, creating a sort of a "traffic jam" , which impedes its motion. The condensed, "traffic jam"- 
like region in the case when the particles density is not explicitly conserved is characterized by 
a stationary density profile which vanishes as an exponential function of the distance from the 
tracer. The characteristic length and the amplitude of the density relaxation function have 
been calculated explicitly. In one dimensional situations with conserved particles density 
the "traffic jam" region grows in size in proportion to the TP mean displacement and no 
stationary profile exists. On the other hand, past the TP the local density is lower than 
the average. In one dimensional situations with conserved particle density this depleted by 
particles region also grows in size in proportion to the TP mean displacement. In higher 
dimensions, we have observed that depending on whether the number of particles is explicitly 
conserved or not, the local density past the TP may tend to the average value at large 
separations from the TP in a completely different fashion: In the non-conserved case the decay 
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of the density is described by an exponential function, while for the conserved particles number 
case it shows an algebraic dependence on the distance, revealing in the latter case especially 
strong memory effects and strong correlations between the particle distribution in the lattice 
gas and the tracer position. Further on, we have found that the terminal velocity V of the TP 
depends explicitly on both the excess density in the "jammed" region in front of the TP, as well 
as on the density in the depleted region past the tracer. We realized that both densities are 
themselves dependent on the magnitude of the tracer velocity, applied external force, as well 
as on the rate of the adsorption/desorption processes and on the rate at which the particles 
can diffuse away of the tracer, which results in effective non-linear coupling between V and E. 
In consequence, in the general case (for arbitrary adsorption/desorption rates and arbitrary 
external force), V can be found only implicitly, as the solution of a transcendental equation 
relating V to the system parameters. This equation simplifies considerably in the limit of a 
vanishingly small external bias; in this case we recover a linear force- velocity relation, akin to 
the so-called Stokes formula. This linear relation signifies that the frictional force exerted on 
the tracer particle by the host medium (the lattice gas) is viscous. The TP mobility, which is 
inverse of the Stokes friction coefficient, thus results from an intricate cooperative behavior. 
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